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1.  Introduction 


Nonlinear  radar  is  well  suited  for  the  detection  of  electronic  devices,  typically  those  containing 
semiconductors  and  whose  traditional  (linear)  radar  cross  sections  are  very  low  owing  to  their 
thin  geometric  profile.  Nonlinear  radar  provides  high  clutter  rejection,1  but  requires  two  major 
tradeoffs:  1)  the  power-on- target  required  to  generate  a  signal-to-noise  ratio  comparable  to  linear 
radar  is  much  higher  than  that  of  linear  radar,2  and  2)  the  radar  must  capture  weak  target 
responses  in  the  presence  of  strong  transmitted  probes.3 

The  nonlinear  radar  studied  in  this  report  transmits  a  single  frequency  at  a  time,  denoted  fo,  and 
receives  twice  this  value,  2/6;  thus,  the  radar  is  harmonic.  Prior  work  has  focused  on  transmitting 
power  at  fo  sufficient  to  generate  a  harmonic  target  response  while  minimizing  the  amount  of 
transmitter-generated  harmonic  that  couples  to  the  receiver.4  More  recent  work  has  focused  on 
amplifying  a  weak  target  response  at  2fo  while  minimizing  the  amount  of  receiver-generated 
harmonic  that  masks  this  response.5  This  report  demonstrates  over-the-air  detection  and  ranging 
of  nonlinear  targets  using  an  ultra- wideband  transmission  and  harmonic  reception. 


2.  Background 


The  basic  architecture  of  a  harmonic  radar  is  shown  in  Fig.  1.  The  radar  transmits  a  single 
frequency  fo,  electromagnetic  energy  at  this  frequency  arrives  at  the  target,  the  target  captures 
and  reradiates  electromagnetic  energy  at  harmonics  of  the  original  frequency,  and  the  radar 
listens  for  a  particular  harmonic,  in  this  case,  2/6.  Reception  of  the  harmonic  indicates  detection 
of  an  electromagnetically  nonlinear  target.  Typically,  because  nonlinear  target  responses  are  so 
weak,  filtering  and  amplification  in  the  receive  chain  are  separated  into  multiple  stages.6-9 


Fig.  1  A  harmonic  radar  that  transmits  at  the  fundamental  frequency  fo  and  receives 
at  the  harmonic  2fo.  Transmitter  =  “Tx”,  receiver  =  “Rx”. 
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In  current  literature,  harmonic  radar  is  typically  narrowband. 1,6-1 5  For  narrowband  systems,  target 
detection  and  classification  are  possible,  but  unless  a  collection  of  antennas  or  a  mobile  platform 
is  implemented,16  range  information  is  not  available.  Ranging  of  nonlinear  targets  is  possible  for 
a  single  stationary  pair  of  antennas,  though,  if  harmonic  responses  are  collected  over  an  ultra¬ 
wide  bandwidth.17 

By  sweeping  the  transmit  frequency  fo  with  known  amplitude  and  phase  across  an  ultra-wide 
bandwidth  and  recording  the  amplitude  and  phase  of  the  harmonic  return  signal  at  each  2/o,  a 
nonlinear  frequency  response  of  the  radar  environment  is  constructed.  An  inverse  Fourier 
transform  of  this  nonlinear  frequency  response  reveals  the  range  to  the  target.  Presented  in  this 
report  is  an  extension  of  a  bench-top,  wireline  experiment  that  we  recently  developed17  to  a 
wireless,  over-the-air  data  collection  using  a  single  ultra-wideband  antenna  for  transmission  and 
reception.  Detection  and  ranging  are  demonstrated  using  3  electromagnetically  nonlinear  targets. 


3.  Experiment 


A  block  diagram  of  the  experimental  ultra-wideband  harmonic  radar  constructed  at  the  Adelphi 
Laboratory  Center  (Building  507)  is  shown  in  Fig.  2.  The  power  budget  for  the  signal  transmitted 
at  fo  and  received  at  2 fo  is  given  in  Fig.  3.  It  should  be  noted  that  this  power  budget  is 
approximate,  because  the  gain/loss  of  each  component/cable  depends  on  fo. 


Fig.  2  Block  diagram  of  an  experimental  ultra-wideband  harmonic  radar 
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Fig.  3  Approximate  power  budget  for  the  experimental  harmonic  radar 

The  signal  source  is  the  Tektronix  AWG7052  arbitrary  waveform  generator  (AWG).  It  outputs 
signals  in  the  frequency  range  fo  =  700  to  900  MHz,  at  a  power  of -3  dBm  (0.5  mW).  The 
remainder  of  the  transmit  chain  is  intended  to  strongly  amplify  fo  and  strongly  attenuate  2fo.  The 
AWG  signal  is  filtered  to  remove  harmonics  using  2  cascaded  MiniCircuits  SLP-1000+  lowpass 
filters.  The  relatively  weak  signal  at  fo  is  increased  by  40  dB  by  the  Amplifier  Research 
AR4W1000  power  amplifier.  The  output  of  the  AR4W1000  is  filtered  by  a  MiniCircuits  NLP- 
1000+  lowpass  filter.  On  its  way  to  the  diplexer  pair,  the  transmitted  signal  is  sampled  using  the 
Hewlett  Packard  1 1607-60001  (20  dB)  directional  coupler.  The  sampled  signal  is  captured  by  1 
channel  of  the  Lecroy  Wavemaster  8300A  oscilloscope.  Most  of  the  transmitted  power  passes 
through  the  coupler  (with  negligible  loss)  to  the  input  of  the  first  diplexer.  Approximately 
33  dBm  (2  W)  reaches  the  antenna  to  be  transmitted  (at  each  frequency). 

Two  diplexers  are  placed  between  the  transceiver  and  the  radar  antenna.  The  lowpass  output  of 
the  first  diplexer  is  connected  to  the  lowpass  output  of  the  second.  The  highpass  output  of  the 
first  diplexer  is  connected  to  a  matched  load  (50  Q),  which  functions  as  a  fourth  filter  in  the 
transmit  path.  Unlike  the  3  MiniCircuits  lowpass  filters,  however,  the  diplexer  terminates 
(absorbs)  harmonics,  i.e.,  it  does  not  reflect  them  back  into  the  transmitter  where  they  may 
couple  into  the  receiver.  The  signal  that  passes  through  the  diplexers  —  in  the  forward  direction 
—  is  sent  to  the  ETS  Lindgren  3164-06  horn  antenna  (“vertical”  port),  and  this  antenna  transmits 
fo  into  the  air.  The  harmonics  generated  by  the  target  are  received  by  the  same  antenna  (on  the 
same  polarization). 

Because  the  highpass  output  of  the  second  diplexer  is  connected  to  the  receive  chain,  only 
harmonics  pass  through  the  diplexers  in  the  reverse  direction.  Therefore,  the  diplexers,  connected 
back-to-back,  1)  provide  lowpass  filtering  in  the  forward  direction,  2)  provide  highpass  filtering 
in  the  reverse  direction,  and  3)  allow  a  single  antenna  to  transmit  fo  and  receive  2 fo.  It  should  be 
emphasized  that  the  Tx  and  Rx  antennas  need  not  be  the  same  antenna.  (If  the  Tx  and  Rx 
antennas  were  separated,  the  lowpass  output  of  the  first  diplexer  would  feed  directly  into  the  Tx 
antenna  instead  of  passing  through  the  second  diplexer,  and  the  lowpass  output  of  the  second 
diplexer  would  be  terminated  in  a  matched  load.  The  transceiver  configuration  would  otherwise 
remain  the  same.) 
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The  receive  chain  (beyond  the  diplexers)  consists  of  4  filters  and  3  amplifiers,  intended  to 
strongly  attenuate  fo,  strongly  amplify  2 fo,  and  ultimately  capture  the  amplitude  and  phase  of 
signals  in  the  frequency  range  2/o  =  1400  to  1800  MHz.  Each  of  the  4  filters  is  a  MiniCircuits 
YHF-1320+  highpass  filter.  The  first  low-noise  amplifier  is  the  MiniCircuits  ZRL-2400LN+;  the 
following  2  amplifiers  are  MiniCircuits  PSA4-5043  evaluation  boards,  powered  by  the  Agilent 
E3631A  (triple-output  DC  power  supply).  The  received  signal  is  captured  by  a  second  channel  of 
the  Lecroy  Wavemaster  8300A  oscilloscope. 

Capturing  the  transmit  and  received  signal  simultaneously  is  necessary  to  determine  the  phase 
difference  between  the  outgoing/)  and  incoming  2fo.  This  phase  difference,  captured  across  a 
wide  band  of  frequencies,  is  critical  to  determining  range-to-target  (as  calculated  in  Section  4). 

S-parameter  measurements  for  portions  of  the  radar  were  recorded  using  the  Agilent  N9923A 
(“FieldFox”)  portable  vector  network  analyzer.  Some  of  the  S-parameters  are  reported  in  this 
section,  but  most  are  reported  in  Appendix  A.  The  magnitude  of  the  transmission  S-parameter 
(S21)  for  both  the  transmitter  (between  the  AWG  and  the  antenna  port)  and  the  receiver  (between 
the  antenna  port  and  the  oscilloscope)  are  given  in  Fig.  4.  The  net  effect  of  the  lowpass  filters, 
power  amplifier,  and  diplexers  (in  the  forward  direction)  is  to  amplify  fo  by  about  36  dB  and 
attenuate  2fo  substantially.  The  net  effect  of  the  diplexers  (in  the  reverse  direction),  highpass 
filters,  and  low-noise  amplifiers  is  to  amplify  2 fo  by  about  50  dB  and  attenuate  fo  substantially. 


400  600  800  1000  1200  1400  1600  1800  2000 

Frequency  (MHz) 


Fig.  4  Overall  measured  l&il  for  the  transmit  and  receive  chains 
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Measured  S-parameters  for  the  diplexer  pair  are  given  in  Fig.  5.  The  input  from  the  transmitter  is 
denoted  Port  1,  the  output  to  the  antenna  is  denoted  Port  2,  and  the  output  to  the  receiver  is 
denoted  Port  3.  Signals  between />  =  700  to  900  MHz  pass  from  Port  1  to  Port  2  with  <2  dB  of 
loss  (|&l|  >  —2  dB  in  this  range),  and  less  than  1 0%  of  the  energy  from  signals  between 
2fo  =  1400  to  1800  MHz  reflect  from  the  diplexer  back  into  the  transmitter  (|5i  i  <  -10  dB  in  this 
range).  Signals  between  2fo  =  1400  to  1800  MHz  pass  from  Port  2  to  Port  3  with  <1  dB  of  loss 
(|&2|  >  -1  dB  in  this  range). 


Fig.  5  Measured  S-parameters  for  the  pair  of  (back-to-back)  diplexers 
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Pictures  of  the  experimental  harmonic  radar  are  shown  in  Figs.  6-1 1 .  Figure  6  is  a  picture  of  the 
signal  generation  and  capture  equipment,  sitting  on  a  cart  and  connected  to  the  Tx/Rx  antenna  by 
a  single  6-ft  SMA  cable.  Figure  7  is  a  closer  view  of  the  Tx/Rx  antenna  and  the  target  location 
(atop  an  empty  cardboard  box,  12  ft  away  from  the  round  metal  back-plane  of  the  antenna). 
Figure  8  is  a  picture  of  the  Panasonic  ToughBook  laptop,  running  Matlab,  used  to  collect  and 
process  the  radar  data.  Figure  9  is  a  closer  view  of  the  directional  coupler  and  the  receive  chain. 
Although  the  diplexers  are  visible  in  Fig.  9,  a  much  closer  view  of  them  is  shown  in  Fig.  10. 
Figure  1 1  is  a  picture  of  the  target  set  —  3  different  nonlinear  radio-frequency  electronic  targets 
(a  connectorized  handheld  radio,  an  unpowered  low-noise  amplifier,  and  a  disconnected  radio¬ 
frequency  mixer)  and  2  different  antennas  (6  inches  long  vs.  16  inches  long). 


Fig.  6  Picture  of  the  experimental  harmonic  radar:  signal  capture  and  generation  cart 
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Fig.  7  Picture  of  the  Tx/Rx  antenna  and  target  emplacement:  Adelphi  Laboratory  Center, 
Building  507,  “sandbox”  area. 


Fig.  8  Picture  of  the  laptop  controller  and  the  power  amplifier 


7 


11607-60001 

coupler 


to/from 

antenna 


iow-noise 

amplifiers 

(x3) 


diplexers 

(x2) 


Fig.  9  Picture  of  the  directional  coupler,  diplexers,  and  receive  chain 


Fig.  10  Close-up  picture  of  the  diplexers:  discrete  components 
are  surface-mounted  on  a  Rogers  43  5 0B  substrate 
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Fig.  11  Picture  of  the  target  set:  MiniCircuits  ZQL- 

1900LNW+  amplifier  (not  powered),  Motorola  FV300 
radio  (antenna  port),  MiniCircuits  ZX05-20H-S+ 
mixer  (intermediate  frequency  [IF]  port) 


4.  Theory  and  Data  Processing 


For  this  experiment,  the  waveform  implemented  to  collect  data  across  the  frequency  range 
/o  =  700  MHz  to  /o  =  900  MHz  was  a  linear  chirp.  If  the  transmitted  waveform  is  a  constant- 
amplitude  sinusoid  whose  frequency  sweeps  linearly  from  /start  to  /  nd,  the  voltage  applied  to  the 
transmit  antenna  may  be  represented  by  the  following  form: 

v—W^.cos^-/.  (<H  =  +  (1) 

1  env 

Tem  is  the  length  of  the  transmitted  pulse  (while  it  is  active).  If  the  radar  receiver  is  designed 
such  that  a  single  harmonic  M  is  selected  for  reception,  the  target  response  —  measured  back  at 
the  radar  —  takes  the  form 

W  =  \VM  I  cos[2;r  •  M  •/„(«)•(«  -  2 r)+  }]  ,,, 
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where  2  r  is  the  time-of- flight  of  the  radar  wave  from  the  transmitter,  to  the  target,  and  back  to 
the  receiver.  By  capturing  the  amplitude  and  phase  of  the  harmonic  response  at  each  frequency, 
the  following  target  response  may  be  constructed: 


f}KL(f) 


V 

rec  _ 

'K  , 

/“/tgt 

V 

trans 

{  0 

/“/tgt 

*BJ2 

>B,J2 


(3) 


//nl  is  the  nonlinear  (harmonic)  frequency  response  of  the  target.  The  response  is  nonzero  over  a 
finite  bandwidth  Btgt  (limited  either  by  the  target  response  or  by  the  transmit  frequency  set), 
centered  at  ftg t.  The  phase  shift  <f>{  Vm]  is  assumed  constant  over  the  band-of-interest  and  is 
therefore  suppressed.  Under  the  assumption  that  |  Vm\  is  also  constant  over  the  target  response 
band,  the  inverse  Fourier  transform  of  T/nl  is 

r  (A  _  S^n{^ tgt  ~  ^r)| I y  \yM-\  y(2^)(/tgt)h-2r) 

NLl  j"  4- 2t)  1  m|  °  f4. 


which  is  a  sine  function,  delayed  by  twice  the  time-of-flight  of  the  received  radar  wave,  and 
scaled  by  the  harmonic  amplitude  response  of  the  target.  Plotted  against  d  =  uPt/2  instead  of  time, 
where  uP  is  the  velocity  of  propagation  of  the  radar  wave,  the  peak  of  /jnl  indicates  the  range  to 
the  nonlinear  target.  For  a  radar  environment  that  is  purely  linear  (i.e.,  when  no  nonlinear  target 
is  present),  /?nl  =  0  for  all  values  of  t. 


5.  Data  Capture  and  Results 


Instrument  control  for  this  experiment  was  performed  using  Matlab  (version  R2013b,  including 
the  Instrument  Control  toolbox,  over  a  general  purpose  interface  bus  [GPIB]-universal  serial  bus 
[USB]  interface).  A  screenshot  of  the  graphical  user  interface  (GUI)  created  to  upload 
waveforms  to  the  Tektronix  waveform  generator  and  download/process  data  from  the  Lecroy 
oscilloscope  is  given  in  Fig.  12.  For  this  experiment,  a  single  waveform  (whose  parameters  are 
illustrated)  was  uploaded  to  the  AWG.  For  each  target  and  distance  away  from  the  Tx/Rx 
antenna,  a  pair  of  data  traces  was  recorded  (vtrans  vs.  ttrans  and  vrec  vs.  tree ,  inside  each  “MAT”  file) 
and  the  harmonic  impulse  response  /?nl  was  computed.  The  resulting  /?nl  traces  are  shown  in 
Figs.  13-16. 
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Fig.  12  GUI  for  harmonic  radar  instrument  control 
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Fig.  13  Data  processed  for  3  targets  with  6-inch  antennas,  placed  6  ft  away  from  the  radar  antenna 
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Fig.  14  Data  processed  for  3  targets  with  16-inch  antennas,  placed  6  ft  away  from  the  radar  antenna 
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Fig.  15  Data  processed  for  3  targets  with  6-inch  antennas,  placed  12  ft  away  from  the  radar  antenna 
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Fig.  16  Data  processed  for  3  targets  with  16- inch  antennas,  placed  12  ft  away  from  the  radar  antenna 

First,  the  6-inch  antenna  was  connected  to  each  target,  and  each  target  was  placed  on  a  cardboard 
box,  6  ft  away  from  the  antenna.  The  antenna  was  pointed  directly  at  the  target  (boresight).  The 
data  recorded  from  the  three  different  targets  (along  with  a  no-target,  “background”  case)  were 
processed  into  Anl  and  appear  in  Fig.  13.  The  16-inch  antenna  was  substituted  for  the  6-inch 
antenna  and  the  same  data  were  recorded  and  processed  into  Fig.  14.  Another  two  data  sets  were 
recorded  with  the  same  3  targets  and  2  antennas  at  a  distance  of  12  ft  away  from  the  antenna  and 
these  data  sets  were  processed  into  Figs.  15  and  16. 

Since  /?nl  is  nearly  zero  for  all  d  in  the  no-target  case,  the  system  noise  is  minimal.  Although  the 
peaks  in  Figs.  13-16  are  somewhat  diffuse  (i.e.,  the  range  resolution,  the  width  of  the  main  lobe 
of  the  response,  is  approximately  5  ft),  the  target  responses  are  easily  (visibly)  discernible  from 
the  system  noise.  The  drop  in  the  peak  response  is  approximately  20  dB  from  d  =  6  ft  to 
d  =  12  ft,  but  the  targets  at  12  ft  are  still  detected.  Detection  and  ranging  of  the  3  nonlinear 
targets  —  using  a  wideband  chirped  waveform,  harmonic  detection,  and  computation  of  a 
nonlinear  impulse  response  —  have  been  demonstrated,  up  to  a  distance  of  12  ft. 
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6.  Conclusion 


An  experimental  harmonic  radar  has  been  constructed  and  successfully  tested  at  transmit 
frequencies  between  700  to  900  MHz  and  an  output  power  of  2  W.  Detection  of  the  nonlinear 
targets  was  accomplished  by  measuring  a  return  signal  (above  the  system  noise)  at  twice  the 
transmit  frequency.  Range-to-target  was  determined  by  sweeping  the  transmit  frequencies  across 
an  ultra-wide  bandwidth,  capturing  the  amplitude  and  phase  of  the  harmonic  return  signal,  and 
computing  the  inverse  Fourier  transform  of  the  wideband  harmonic  response.  Detection  and 
ranging  of  3  nonlinear  targets  (radio-frequency  electronic  components  connected  to  straight- wire 
antennas)  were  demonstrated,  up  to  a  distance  of  12  ft  away  from  the  radar  antenna. 
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Appendix  A.  Measured  S-Parameters  for  Each  Component  of 

the  Harmonic  Radar 


Figures  A-l  through  A-8  show  the  measured  S-parameters  for  each  component  of  the  harmonic 
radar . 


Fig.  A-l  Measured  S-parameters  for  the  MiniCircuits  SLP-1000+  lowpass  filter  pair 


Fig.  A-2  Measured  S-parameters  for  the  Amplifier  Research  AR4W1000  power  amplifier 
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Fig.  A-3  Measured  S-parameters  for  the  MiniCircuits  NLP-1000+  lowpass  filter 


Fig.  A-4  Measured  S-parameters  for  the  HP  11607-60001  directional  coupler: 
port  1  =  input,  port  2  =  through  (output),  port  3  =  coupled. 
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Fig.  A-5  Measured  S n  for  the  ETS  Lindgren  3164-06  quad-ridge  horn  antenna 


Fig.  A-6  Measured  S-parameters  for  the  MiniCircuits  VHF-1320+  highpass  filter 
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Fig.  A-7  Measured  S-parameters  for  the  MiniCircuits  ZRL-2400LN+  low-noise  amplifier 


Fig.  A-8  Measured  S-parameters  for  the  MiniCircuits  PSA4-5043  low-noise  amplifier 
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Intentionally  Left  Blank. 
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Appendix  B.  Matlab  Script  to  Control  the  Graphical  User  Interface 


23 


o  o  o 
o  o  o 


function  radar_nonlinear_OpeningFcn (hOb j ect ,  eventdata,  handles,  varargin) 


%  connect  to  the  AWG  (Tektronix  AWG7052) 

%  and  to  the  oscilloscope  (Lecroy  Wavemaster  8300A) 
global  awg  osc; 

awg  =  visa (  '  ni  '  ,  ' GPIBO :  : 2 :  : INSTR ' )  ; 
osc  =  visa (  ' ni ' , ' GPIBO : : 6 : : INSTR ' ) ; 

f open  (awg) 

fwrite (awg,  ’ AWGCONTROL : D0UTPUT1 : STATE  1’);  %  turn  'direct  output'  on 

fclose (awg) 

%  define  strings  for  the  Tx  channel  (C2),  Rx  channel  ( C 3 ) , 

%  Tx  channel  averaged  (F2),  and  Rx  channel  averaged  (F3) 

global  data_trc_l  data_trc_2  avg_trc_l  avg_trc_2  sweeps; 

data_trc_l  =  'C3'; 

data_trc_2  =  'C4'; 

avg_trc_l  =  ' F3 ' ; 

avg_trc_2  =  '  F4  '  ; 

sweeps  =  1; 


%  turn  C3  and  C4  off;  turn  F3  and  F4  on 


f open (osc) 

fprintf (osc, [avg  trc  1  ' 
num2str (sweeps) ] ) ; 

: DEF  EQN, ' 'AVG ( ' 

data  trc  1 

fprintf (osc, [avg  trc  2  ' 
num2str (sweeps) ] ) ; 

: DEF  EQN, ' 'AVG ( ' 

data  trc  2 

fprintf (osc, [data  trc  1 

' : TRACE  OFF ' ] ) ; 

fprintf (osc, [data  trc  2 

' : TRACE  OFF ' ] ) ; 

fprintf (osc, [avg  trc  1  ' 

: TRACE  ON ' ] ) ; 

fprintf (osc, [avg  trc  2  ' 

: TRACE  ON ' ] ) ; 

fprintf (osc, [avg  trc  1  ' 

: FRST ' ] ) ; 

%  reset 

fprintf (osc, [avg  trc  2  ' 
fclose (osc) 

: FRST ' ] ) ; 

' ) , SUMMED' ' , SWEEPS, ' 
' ) , SUMMED' ' , SWEEPS, ' 


the  averaged  sweeps 


handles . output  =  hObject; 
guidata (hOb j ect ,  handles); 


ooooooooooooooooooooooooooooooooooooo 

%%%  generate  single-tone  RF  pulse  %%% 

ooooooooooooooooooooooooooooooooooooo 


function  pushbuttonl_Callback (hOb j ect ,  eventdata,  handles) 
global  awg; 

%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . edit6, ' String ') ; 
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%  user-defined  pulse  power  while  the  pulse  is  active 
global  P_env; 

P_env  =  str2num (get (handles . editl ,  ' String ')) ; 

%  user-defined  carrier  frequency  for  the  RF  pulse 
global  f_pulse; 

f_pulse  =  str2num (get (handles . edit2 ,  ’ String ')) ; 

%  user-defined  time  for  the  RF  pulse  to  be  active 
global  T_env; 

T_env  =  str2num (get (handles . edit4 ,  ' String ')) ; 

%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  =  str2num (get (handles . edit5 ,' String ')) /100  ; 

%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 ,  ’ String ')) ; 

%  generate  the  single-tone  pulse  waveform 
%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  fsample; 

[v, f_sample]  =  awg_pulse ( f_pulse, T_envA P_env, duty_cyclef  f_sample_max) 
%  awg_plot (v, f_sample) 

oooooooooooooooooooooooooooooooooooooooooooooooooo 


global  samples; 
samples  =  length (v) ; 

%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros (1, samples) ; 

marker_l ( 1 : round ( samples/10 ) )  =  1; 

marker_2  =  zeros (1, samples) ; 

marker_2 ( 1 : round ( samples/10 ) )  =  1; 

global  buffer; 
buffer  =  samples; 

set (awg, ' InputBuf f erSize ' , buffer) ; 
set (awg, ' OutputBuf f erSize ’ , buffer) ; 

%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 

%  set  the  AWG  sampling  frequency,  select  the  waveform  to  be  played, 
%  and  turn  the  AWG  output  on 
f open (awg) ; 

fwrite (awg,  'AWGCrDOUTl  ON ' ) ; 

fwrite (awg,  [* S0URCE1 : FREQUENCY  '  num2 str ( f_sample/ 1 0 A 9 )  '  GHZ']); 
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fwrite (awg,  [' SOURCE 1 : WAVEFORM  wf_name  ' 

fwrite (awg,  'OUTPUT  ON'); 
fwrite (awg,  ' AWGCONTROL : RUN ' ) ; 
fclose (awg) ; 


oooooooooooooooooooooooooooooooooooooo 


%%%  generate  linear  RF  chirp  pulse  %%% 


function  pushbutton2_Callback (hOb j ect ,  eventdata,  handles) 
global  awg; 


%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . editl2 String ') ; 

%  user-defined  chirp  power  while  the  pulse  is  active 
global  P_env; 

P_env  =  str2num (get (handles . edit8 ,  ' String ')) ; 

%  user-defined  start  frequency  for  the  chirp 
global  f_start; 

f_start  =  str2num (get (handles . edit 9 String ')) ; 

%  user-defined  end  frequency  for  the  chirp 
global  f_end; 

f_end  =  str2num (get (handles . editl3 r  ’ String ')) ; 

%  user-defined  length  of  the  chirp 
global  T_env; 

T_env  =  str2num (get (handles . editlO ,  ’ String ’)) ; 

%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  =  str2num (get (handles . editl 1 , ' String ')) /100 ; 

%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 String ')) ; 

%  generate  the  RF  chirp  waveform 

%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  f_sample; 

[v, f_sample]  =  awg_chirp ( f_start ,  f_end^  T_envf  P_envf  duty_cycleA  f_sample_max) 
%  awg_plot (v, f_sample) 


global  samples; 
samples  =  length (v) ; 
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%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros ( 1 , samples) ; 

marker_l ( 1 : round ( samples/2 0 ) )  =  1; 

marker_2  =  zeros (1, samples) ; 

marker_2 ( 1 : round ( samples/2 0 ) )  =  1; 

global  buffer; 
buffer  =  samples; 

set (awg, ' InputBuf f erSize ' , buffer) ; 
set (awg, ' OutputBuff erSize ’ , buffer) ; 

%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 

%  set  the  AWG  sampling  frequency,  select  the  waveform  to  be  played, 
%  and  turn  the  AWG  output  on 
f open (awg) ; 

fwrite (awg,  'AWGCrDOUTl  ON ' ) ; 

fwrite (awg,  [' S0URCE1 : FREQUENCY  '  num2 str ( f_sample/ 1 0 A 9 )  '  GHZ ' ] ) ; 

fwrite (awg,  [’ SOURCE 1 : WAVEFORM  " '  wf_name  ""]); 
fwrite (awg,  'OUTPUT  ON'); 
fwrite (awg,  ' AWGCONTROL : RUN ' ) ; 
f close (awg) ; 


%%%  generate  multitone  RF  pulse  %%% 

ooooooooooooooooooooooooooooooooooo 


function  pushbutton4_Callback (hObj ect ,  eventdata,  handles) 


global  awg; 


%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . editl8 ,' String ') ; 
global  N; 

N  =  str2num (get (handles . editl 9 , 'String')); 


%  user-defined  power  per  tone  while  the  pulse  is  active 
global  P_tone; 

P_tone  =  str2num (get (handles . editl4 ,' String ')) ; 

%  user-defined  center  frequency  for  the  tones 
global  f_c; 

f_c  =  str2num (get (handles . editl5 ,' String ')) ; 

%  user-defined  time  for  the  multitone  pulse  to  be  active 
global  T_env; 

T_env  =  str2num (get (handles . editl 6 ,' String ')) ; 
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%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  =  str2num (get (handles . editl7 ,  ' String ')) /100 ; 

%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 ,  ' String ')) ; 

%  generate  the  multitone  waveform 

%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  fsample; 

[v ,  f_sample]  =  awg_multitone (N, f_c,  T_env,  P_tone, duty_cycle, f_sample_max) 
%  awg_plot (v, f_sample) 


global  samples; 
samples  =  length (v) ; 

%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros (1, samples) ; 

marker_l ( 1 : round ( samples/10 ) )  =  1; 

marker_2  =  zeros (1, samples) ; 

marker_2 ( 1 : round ( samples/10 ) )  =  1; 

global  buffer; 
buffer  =  samples; 

set (awg, ’ InputBuf f erSize ’ , buffer) ; 
set (awg, ' OutputBuf f erSize ' , buffer) ; 

%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 

%  set  the  AWG  sampling  frequency,  select  the  waveform  to  be  played, 
%  and  turn  the  AWG  output  on 
f open (awg) ; 

fwrite (awg,  'AWGCrDOUTl  ON'); 

fwrite (awg,  [’ SOURCE 1 : FREQUENCY  '  num2 str ( f_sample/ 1 0 A 9 )  '  GHZ ' ] ) ; 

fwrite (awg,  [’ SOURCE 1 : WAVEFORM  wf_name  '"']); 

fwrite (awg,  'OUTPUT  ON'); 
fwrite (awg,  ' AWGCONTROL : RUN ' ) ; 
f close (awg) ; 


oooooooooooooooooooooooooooooooooooooooo 


%%%  generate  stepped-f requency  pulse  %%% 


function  pushbutton5_Callback (hOb j ect ,  eventdata,  handles) 
global  awg; 
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%  user-defined  waveform  name  to  be  uploaded  to  the  AWG 
global  wf_name; 

wf_name  =  get (handles . edit24 ,  ' String '  ) ; 

%  user-defined  power  at  each  frequency  while  the  pulse  is  active 
global  P_env; 

P_env  =  str2num (get (handles . edit2 0 String ')) ; 

%  user-defined  start  frequency  for  the  steps 
global  f_start; 

f_start  =  str2num (get (handles . edit2 6 ,  1 String '))  ; 

%  user-defined  end  frequency  for  the  steps 
global  f_end; 

f_end  =  str2num (get (handles . edit2 1 ,  ’ String '))  ; 

%  user-defined  frequency  step  size 
global  f_delta; 

f_delta  =  str2num (get (handles . edit25 ,  ' String '))  ; 

%  user-defined  time  for  the  stepped  waveform  to  be  active 
global  T_env; 

T_env  =  str2num (get (handles . edit22 ,  ’ String ')) ; 

%  user-defined  duty  cycle 
global  duty_cycle; 

duty_cycle  =  str2num (get (handles . edit23 String ')) /100 ; 

%  user-defined  maximum  sampling  rate  for  the  AWG 
global  f_sample_max; 

f_sample_max  =  str2num (get (handles . edit3 ,  ’String')); 


%  generate  the  stepped-f requency  waveform 
%  and  compute  an  appropriate  AWG  sampling  frequency 
global  v  f _s ample; 

[v,f_sample]  = 

awg_stepped ( f_start ,  f_end, f_delta, T_envf  P_envf  duty _ cycle ^ f_sample_max) 

%  awg_plot (v, f_sample) 


global  samples; 
samples  =  length (v) ; 

%  provide  two  markers  (triggers)  in  sync  with  the  AWG  output 

marker_l  =  zeros (1, samples) ; 

marker_l ( 1 : round ( samples/2 ) )  =  1; 

marker_2  =  zeros (1, samples) ; 

marker_2 ( 1 : round ( samples/2 ) )  =  1; 

global  buffer; 
buffer  =  samples; 
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set (awg, 1 InputBuf f erSize ' , buffer) ; 
set (awg, 1 OutputBuff erSize 1 , buffer) ; 


%  convert  the  real  waveform  data  to  binary 
%  and  upload  the  binary  data  to  the  AWG  for  playback 
[binblock,  header,  bytes]  =  awg_binary (v,  marker_l,  marker_2); 
awg_upload (awg,  buffer,  wf_name,  binblock,  header,  bytes); 


%  set  the  AWG  sampling  frequency, 
%  and  turn  the  AWG  output  on 
f open (awg) ; 

fwrite (awg,  'AWGCrDOUTl  ON'); 
fwrite (awg,  [ 1 S0URCE1 : FREQUENCY  1 
fwrite (awg,  [' SOURCE 1 : WAVEFORM  "' 
fwrite (awg,  'OUTPUT  ON'); 
fwrite (awg,  ' AWGCONTROL : RUN ' ) ; 
f close (awg) ; 


select  the  waveform  to  be  played. 


num2str (f_sample/10A9)  '  GHZ']); 
wf  name  '"']); 


%%%  clear  waveforms  &  reset  AWG  %%% 


ooooooooooooooooooooooooooooooooooo 


function  pushbutton3_Callback (hOb j ect ,  eventdata,  handles) 

global  awg; 

f open (awg) ; 
fwrite (awg, ' *RST ' ) ; 
fwrite (awg, ' *CLS ' ) ; 
pause ( 1 ) 

fwrite (awg,  ' AWGCONTROL : D0UTPUT1 : STATE  1'); 
f close (awg) ; 


ooooooooooooooooooooooooooooooooooooooo 


%%%  capture  and  process  Rx  waveform  %%% 


function  pushbutton7_Callback (hOb j ect ,  eventdata,  handles) 
global  awg  osc; 
global  trig_src; 

switch  get (handles . popupmenu2 , 'Value') 
case  1 


trig  src  = 
case  2 

'  EX '  ; 

o 

o 

trigger  = 

external 

trig  src  = 
case  3 

'Cl '  ; 

o 

o 

trigger  = 

channel  Cl 

trig  src  = 
otherwise 

trig  src  = 

'  C2  '  ; 

'EX'; 

o 

o 

trigger  = 

channel  C4 

end 
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%  user-defined  trigger  level 
global  trig_lev; 

trig_lev  =  str2num (get (handles . edit37 ,  ' String ')) *10 A-3 ; 

%  user-defined  total  time  for  Tx  and  Rx  data  to  be  collected 
global  t_total; 

t_total  =  str2num (get (handles . edit32 ,  ' String ')) *10A-6; 
t  div  =  t  total/10; 


%  user-defined  voltage  per  division  visible  on  the  oscilloscope 
global  v_div_l  v_div_2 ; 

v_div_l  =  str2num (get (handles . edit33 ,' String '))  *10 A-3 ; 
v_div_2  =  str2num (get (handles . edit34 ,  '  String ')) *10 A-3 ; 

%  user-defined  voltage  offset  on  the  oscilloscope 
global  v_offset; 
v_offset  =  Oe-3; 

%  set  trigger  delay  to  1/2  of  the  data  collection  time  so  that 
%  the  maximum  amount  of  data  is  collected  AFTER  the  trigger  point 
global  trig_delay; 
trig_delay  =  -t_total/2; 

%  user-defined  Matlab  file  to  store  Tx  and  Rx  time  &  voltage  vectors 
global  data_name; 

data_name  =  get (handles . edit36, ' String ') ; 
global  trig_mode; 

trig_mode  =  'SINGLE';  %  use  a  single  trigger  for  each  sweep 

global  f_capture; 
f_capture  =  10e9; 


global  buffer; 

buffer  =  10*t_div*f_capture  +  26; 
set (os c,  ' InputBuf f erSize ' , buffer) ; 
set (os c,  ' OutputBuff erSize ' , buffer)  ; 

global  data_trc_l  data_trc_2  avg_trc_l  avg_trc_2  sweeps; 
sweeps  =  str2num (get (handles . edit38 ,  'String')); 


set  the  oscilloscope  to  record  data  according  to  user-input  values 


f open (osc) ; 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc, 
fprintf (osc. 


[ ' TRSE  EDGE, SR,  '  trig_src  '  ,  HT,  OFF ' ] )  ; 

[ ' TRMD  '  trig_mode] ) ; 

[ ' TRDL  '  num2str ( trig_delay/10 A-6)  '  US']); 

[ ' EX : TRLV  '  num2 str ( trig_lev/ 1 0 A-3 )  '  MV']) 

[ ' TDIV  '  num2str(t  div/10A-6)  '  US ' ] ) ; 


[data_trc_l 
[data_trc_2 
[data  trc  1 


VDIV  '  num2str (v_div_l/10 A-3 ) 
VDIV  '  num2str (v_div_2/10A-3) 
OFST  '  num2str(v  of f set/10 A-3 ) 


MV ' ] ) ; 
MV ' ] ) ; 

'  MV ' ] )  , 
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MV’ ] ) ; 


fprintf (osc,  [data_trc_2  ':OFST  '  num2str (v_off set/10 A-3 )  ' 

f close (osc) / 

%  turn  the  appropriate  data  traces  on/off 
%  and  set  two  of  the  traces  to  each  average  the  Tx  or  Rx  data  stream 
f open  (osc) ; 

fprintf (osc, [data_trc_l  ’ : TRACE  OFF' ] ) ; 
fprintf (osc, [data_trc_2  ' : TRACE  OFF' ] ) ; 


fprintf (osc, [avg  trc  1  ' 
fprintf (osc, [avg  trc  2  ' 
fprintf (osc, [avg  trc  1  ' 
fprintf (osc, [avg  trc  2  ' 

TRACE  ON' 

TRACE  ON' 
FRST ' ] ) ; 
FRST ' ] ) ; 

] )  ; 

] )  ; 

fprintf (osc, [avg  trc  1  ' 
num2str (sweeps) ] ) ; 

DEF  EQN, ' 

' AVG ( ' 

data  trc  1 

' ) , SUMMED' ' , SWEEPS, 

fprintf (osc,  [avg  trc  2  '  : 
num2str (sweeps) ] ) ; 

: DEF  EQN, ' 

' AVG ( ' 

data  trc  2 

' ) , SUMMED' ' , SWEEPS, 

for  counter  =  1 : sweeps 
fprintf (osc,  ' *TRG '  )  ; 
fprintf (osc, ' WAIT ' ) ; 

end 

f close (osc) ; 

%  download  the  oscilloscope  data  into  Matlab 

[v_trans,  v_rec,  t_trans,  t_rec]  =  osc_capture  (  osc,  buffer,  avg_trc_l, 
avg_trc_2  ) ; 

%  stop  the  waveform  generator  playback 
f open (awg) ; 

fwrite (awg,  'OUTPUT  OFF'); 
fwrite (awg,  ' AWGCONTROL : STOP ' ) ; 
f close (awg) ; 

%  save  the  raw  Tx  and  Rx  vectors  to  a  MAT  file 
save (data_name, ' t_trans ' , ' v_trans ' , ' t_rec ' , ' v_rec ' ) ; 

%  process  and  display  the  Tx  and  Rx  data 
osc_process (  v_trans,  v_rec,  t_trans,  t_rec  ) ; 


32 


Appendix  C.  Matlab  Function  to  Generate  a  Chirp  Waveform 
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function  [v,  f_sample]  =  awg_chirp (  f_start,  f_end,  T_env,  P_env,  duty_cycle, 
f_sample_max  ) 

f_res  =  1/T_env;  %  frequency  resolution  (Hz) 


f_sample  =  2*f_res; 

while  (f_sample  <=  f_sample_max/2 )  %  adjust  AWG  sample  rate 

f_sample  =  2*f_sample; 

end 

f_start  =  round ( f_start/f_res ) *f_res ;  %  adjust  start  frequency 

f_end  =  round ( f_end/f_res ) *f_res ;  %  adjust  end  frequency 

oooooooooooooooooooooooooooooooooooooooooooooooooo 


t  sbb  =  1/f  sample; 

o 

o 

sampling  time 

t  =  0:t  sbb : T  env-t  sbb; 

o 

o 

time  vector 

A  =  sqrt (10A (P_env/10) *2*50*10A-3) ; 

o 

o 

amplitude  of  chirp  when  ON 

k  =  (f  end-f  start) /T  env; 

o 

o 

chirp  rate  (Hz/s) 

x  =  A  *  cos  (2*pi* (f  startf (k/2 ) *t )  . *t ) ; 

o 

o 

chirp  equation 

if  (duty  cycle  <=  0)  | |  (duty  cycle  >=  1) 

v  =  x  ; 

o 

o 

continuous  wave 

else 

v  =  [x  zeros ( 1 ,  (round ( l/duty_cycle) -1 ) ^length (x) )] ;  %  pulsed 

end 
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Appendix  D.  Matlab  Function  to  Generate  a  Step-Frequency  Waveform 
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function  [v,  f_sample]  =  awg_stepped(  f_start,  f_end,  f_delta,  T_env,  P_env, 
duty_cycle,  f_sample_max  ) 

f_res  =  1/T_env;  %  frequency  resolution  (Hz) 

f_sample  =  2*f_res; 

while  (f_sample  <=  f_sample_max/2 )  %  adjust  AWG  sample  rate 

f_sample  =  2*f_sample; 

end 

freq  =  ( f_start : f_delta : f_end-f_delta) ;  %  instantanous  frequencies 

N  =  length ( freq) ; 

%  adjust  instantaneous  frequencies 
for  counter  =  1 : length ( freq) 

freq ( counter )  =  floor (freq (counter) /f_res) *f_res; 

end 


t_sbb  =  l/f_s ample; 
t  =  0 : t_sbb : T_env-t_sbb; 

A  =  sqrt (10A (P_env/10) *2*50*10A-3)  ; 


%  sampling  time 
%  time  vector 

%  amplitude  of  chirp  when  ON 


index  =  ceil ( (1 : length (t) ) / (length (t) /N) ) ; 
f  =  freq(index); 
x  =  A  *  cos  (2*pi* (f . *t) ) ; 


if  (duty_cycle  <=  0)  | |  (duty_cycle  >=  1) 

v  =  x;  %  continuous  wave 

else 

v  =  [x  zeros ( 1 ,  (round ( l/duty_cycle) -1 ) ^length (x) )] ;  %  pulsed 

end 
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Appendix  E.  Matlab  Function  to  Process  the  Received  Waveform  into  /inl 
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function  osc_process (  v_l,  v_2,  t_l,  t_2  ) 

v_in  =  v_l ; 
v_ref  =  v_in.A2; 
v_out  =  v_2 ; 

t  =  ( 0 : 1 : length (v_in) -1 )  *  le-10; 

v_out  =  bandpass (  t,  v_out,  700e6*2,  900e6*2  ); 

delta_t  =  t(2)-t(l); 
f_sample  =  l/delta_t; 

T  =  t  (length (t) )  +  l/f_sample; 

freq  =  -f_sample/2 : 1/T : ( f_sample/2 ) -1/T; 

V_in  =  f ftshift (f ft (2* (v_in) ) /length (v_in) ) ; 

V_ref  =  f ftshift (f ft (2* (v_ref) ) /length (v_ref) ) ; 

V_out  =  f ftshift (f ft (2* (v_out) ) /length (v_out) ) ; 

H  =  V_out . /V_ref ; 

h  =  ifft  (fftshift (H/2) *length(H) ) ; 

%  adjust  the  velocity  of  propagation  by  dielectric  constant 
dielectric_PTFE  =  2.1; 

velocity_PTFE  =  0.9836e9  /  sqrt (dielectric_PTFE) ; 
d  =  velocity_PTFE  *  t/2; 

oooooooooooooooooooooooooooooooooooooooooooooooooo 


figure  ( 10 ) 

plot (d, real (h) , 1 kf ) 

axis ([0400  Inf] ) 

ylabel ( ' { \ith }_{NL }  (rescaled)') 

xlabel (' Distance  to  Target  (ft)  ') 
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